Preprint INR-889/95 



COMPUTATION OF COEFFICIENT FUNCTIONS FOR AN EXPANSION 
OF PARAPOSITRONIUM DECAY WIDTH USING REDUCE * 



A.A.PIVOVAROV 



Institute for Nuclear Research 
Moscow, 117312, Russia) 



A kind of asymptotic expansion for calculation of the (para)positronium decay rate is 
formulated. The method is completely relativistic and gauge invariant. The expansion 
goes over the inverse square of relative momentum k of two photons in Euclidean region 
and the actual width is determined at k 2 = —m? . m being the electron mass. 



* Talk presented at Fourth AIHENP95 Workshop, Pisa, Italy, 3-8 April 1995 
^E-mail: aapiv@ms2.inr.ac.ru 



2 Computation of coefficient functions 



At present, computation of physical observables with high precision is of great im- 
portance for testing theory. In some cases the accuracy of experimental data is 
better than theoretical predictions that requires further calculations for adequate 
comparison. The standard perturbation theory can then become fairly difficult to 
live with. Thus, the most recent data on lifetime of (ortho)positronium 1,2 seem to 
contradict the present theoretical calculations 3-5 done in the ncxt-to-lcading order 
in fine structure constant a. Full two-loop calculations are not available yet and 
the discrepancy is often parametrized by the value of second order coefficient (a 
factor in front of (a/ir) 2 ) that should be about 250 (e.g. ). This number is consid- 
ered quite big though not inconceivable especially bearing in mind the large first 
coefficient. Despite all speculations one cannot claim the real inconsistency till the 
reliable computation of the second order coefficient and confirmation of experimen- 
tal measurements. Recently another special set of diagrams has been analytically 
computed in 6 . New formulation based on explicit use of nonrelativistic effective the- 
ory has been also implemented for getting a representation that stresses the splitting 
of momentum region into relativistic and nonrelativistic parts 7 ' 8 . Along with these 
efforts to have new additional terms computed there are more formal attempts to 
resum the known series using different kinds of summation technique, for instance, 
Pade approximation. Indeed, rewriting the scries in equivalent form up to relative 
order a one can make the required higher order coefficients essentially smaller 9 . A 
remark is in order here. Strictly speaking the rate should note be represented as a 
power series in a - even some analytical functions of a connected with exact treat- 
ment of an initial approximation can be kept in their concise form. This will not 
mean any violation of consistency because a is not an expansion parameter. 10 For 
fitting experiment, it implies that the change of zero order approximation in Bcthc- 
Salpcter approach could bring new terms of any order in a that can stand on the 
legal ground. Still, it is difficult to understand the required large value of the second 
order coefficient. The technique based on the Bethe-Salpeter equation was widely 
used for computation of quarkonium decay rates 11 ' 12 where a new nonrelativistic 
theory based approach has been recently developed as well. 8,13 

The present note describes a technique of calculating the parapositronium decay 
rate into photons that is inspired by methods of physics of strong interaction. 14 
Results of explicit calculation of parapositronium width up to the first order in the 
fine structure constant are reported. 15 

The Lehmann-Symanzik-Zimmermann reduction 16 leads to the following expres- 
sion for the two-photon decay rate of positronium 

(l(kih(k 2 )\P(p)) = i{2^f5{p -h- k 2 )e 2 e"(k 1 )e"(k 2 )(0\Il^(k)\P(p)), 
T^{x) = iTJ^x/2)J v (-x/2), r M „(fc) = i J TJ^x/2)J v (-x/2)e lkx dx (1) 

where \P(p)) is an approximation for the positronium bound state, p 2 = m 2 Pl k — 
(ki — k2)/2 1 J M is the fermion electromagnetic current, e M (fci) is a polarization 
vector of a photon. Because of momentum conservation, p = k\ + k 2l and mass- 
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shell conditions for the photons, fcf = fcf = 0, there exist kinematical constraints 
on the variables k and p, kp — and k 2 = —mp/i, mp being the "positronium" 
mass. 

The main thing to analyze is therefore a T-product of two electromagnetic cur- 
rents T^(k) that, being sandwiched between the vacuum and positronium state, 
gives a matrix element of the decay process 

(0|T^(A)|P(p)> = e^fcVi^ 2 ), r(P - 77 ) = |T(-^)| 2 ^. (2) 

Little can be said about the vector \P(p)) besides that its quantum numbers are 
J PC = ^. It is quite analogous to the neutral pion state in QCD. To the leading 
order in a Eq. (||) gives its decay width but corrections cannot be unambiguously 
computed unless the electromagnetic contribution to the pion wave function itself is 
properly taken into account. As for the T-product of currents, it can be expanded 
in x when x goes to zero. 17 For the matrix element, Eq. (jl|), this expansion turns 
into a series in the ratio of two scales connected with soft and hard parts of the 
process respectively. The first one is the Bohr radius of the positronium ground state 
rg 1 /am that determines an effective size of the bound state wave function and 
the second is the cutoff length set by the fermion mass m, r cut ~ 1/m. Eventually, 
it is an expansion in r cut /rs, or in a, but not completely since each term can 
contain some higher powers of a as well. The decay rate given by the form factor 
F(k 2 ) will be parametrized with some numbers that describe the inner structure of 
positronium and can not be found without analyzing its dynamics as a bound state. 

To obtain the decay rate one needs the form factor F(k 2 ) only at the physical 
point k 2 = —rap/A. The asymptotic expansion of the T-product can be done for any 
value of k 2 in Euclidean region k 2 < because one does not encounter any physical 
singularities. Fortunately, the physical point belongs to the allowed region and one 
can compute the form factor at the desired point k 2 = —m 2 p/A. The convergence of 
the expansion, which is asymptotic one, at this particular point is not guaranteed. 

Generating formula for the expansion reads 

T^{x) = ^(x/2) 7m G(x) 7 ^(-x/2) + ^{-x/2) lv G{-x)^{x/2) + ... (3) 

where G(x) stands for a full fermion propagator. Since one needs only operators 
with parapositronium quantum numbers the T-product can be antisymmetrized 
with respect to indices fi and v (see Eq. (0)). This is implicitly implied in the 
following. 

The leading term has the form 

T o = «e^y (4) 
and gives the contribution to the form factor F(k 2 ) 



(•5) 
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where fp is defined by the relation (O|'07 Al 75'0|-P(p)) = ip^fp, and a notation k 2 = 
—k 2 = mp/A is used. Clearly, this quantity is similar to the pion decay constant. 
After inserting a correction to the fermion propagator 

AG{k) = (m 2 -fc 2 ) 2 + 2(m 2 -fc 2 ) 2 ' F *»-2**»* F 

into Eq. (Q) and set arguments of fermion fields equal to zero the following result 
for a dimension-five operator connected with a photon emission follows 



TL{k) = 



i , _ _ 2ie fl TV xk p eipF frr ^ x ip 
(m 2 + k 2 ) 2 



Defining the S 2 parameter by the relation 

(0\e^F" T ^\P(p)) - (P P 9 rX -p T 9 pX )ifp5 2 /3 

one finds 

F F( k 2 ) = - 

V ' 3(m 2 + n 2 ) 2 ' 

The quantity 5 2 is obviously of order a at least. 

Next terms of the series are given by derivatives of fermion fields after substitu- 
tion of a free propagator S(x) for the full one G(x). The first-order derivative gives 
a contribution 

j (m 2 + K 2)2 

that turns into 

3 (m 2 + k 2 ) 2 

after averaging between the vacuum and positronium states. This term is also 
suppressed by a due to presence of the difference 4m 2 — m 2 P . 

The second-order derivative leads to an operator ■07 A 75_D Q D i a'0, the matrix 
element of which can be expressed through three numbers A, B, C 

(0$B°S Vts^OO) = *fpAg al3 P x + B(g aX P f3 + / V) + CpW 
though only one of them contributes 

fpA 



F 2 {k z ) = 



(TO 2 + K 2 ) 2 TO 2 + K 2 



I failed to show in general that the quantity A is of the order a in comparison with 
the leading term, Eq. (5). 

The final representation up to dimension-five operators reads 

ptf) 2 -fr g 62 4m 2 -to 2 . A 3m 2 + K 2 

V ' TO 2 + K 2 3(TO 2 + K 2 ) 6(TO 2 + K 2 ) 2(TO 2 + K 2 ) TO 2 + K 2 ' V ' 



Computation of coefficient functions 5 



For a consistent analysis perturbative corrections to the coefficient function of 
the leading operator, Eq. (0), have to be found. Being written in the form 



^ (ib ) = w^y^ a-H"^)), (7) 



they are 



, n Ci . o 12 3 . Ct , o 12 . 

C(«) = + 10z + 2 + 7 + -) + -W( z - 7z ~ 32 + — ) 

3 2 + 1 Z 3 z 

+ 1 -(-2zi + 19z- 1 -± + ^-+*) (8) 
3 2 z + 1 z 



where z = k 2 /m 2 and 



(z + l)ln(z + 1) 



C\ = I dx ]n(zx(l — x) + x) = 2, 

Jo z 

c — Z f 1 ^ X — ^ In a + ^ a — a/ 1 + 4/z (9) 

2 7 2a;(l — x) + 1 a a— 1' 

The electron mass is the pole mass. The actual calculation was performed in the 
M 5-scheme of renormalization with its natural mass parameter mjjg. A relation 
between these two up to one-loop order is 18 

(a fj 2 \ 
1 H (3 m — 5 h 1) . 
W m MS J 

The quantity c(z) does not depend on the renormalization scheme if the physical 
pole mass is used. The program of symbolic manipulations REDUCE has been 
heavily used for computation. 

To verify Eq. (8) one can investigate two simple limits. First, z — > oo that 
corresponds to m — > or k 2 — > oo at fixed m. In this limit the logarithmic z- 
dependence of function c(z) must completely drop out because of the zero anomalous 
dimension of the coefficient function. The result is c(oo) = —7. This result has been 
also checked by independent calculation with taking the massless limit from the 
very beginning. Though contributions of different diagrams have been considerably 
rearranged a full agreement for the final answer was found. The opposite limit 
z — > also must not reveal any singularity. It is easy to see that Eq. (8) is valid, 
c(0) = 7/2. 

Neglecting the difference between 4m 2 and m|>, which is of order a 2 , one gets 

a , . a 
1 + — c 1 =1 + 1.012—. 

An W 4-7T 

Now we briefly consider parameters of Eq. (6). The amplitude fp is connected 
with that of the positronium wave function at the origin 19 . The 8 2 parameter is 
reduced to the matrix element (O|e - 0F M „7 iy ?/'|P(j))) = ip^fpS 2 that is closely related 
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to the corresponding S 2 parameter in strong interaction 20 . As for the A parameter, 
I did not succeed in expressing it in any transparent manner. Both S 2 and A 
parameters related to the positronium wave functions of nonleading twists. 

Eq. (6) represents the main result. One sees that the two-photon decay width of 
parapositronium is expressed through three parameters fp, S 2 and A that require 
a special model to determine them. Unambiguous physical meaning (up to a cor- 
responding a order within perturbative QED) can be given to physical quantities 
that do not contain these parameters. They could be computed independently if 
one knows the positronium wave function from other sources (for instance, lattice 
calculation). The main problem how to define properly the vector \P(p)) cannot be 
answered within the approach; so, strictly speaking, the model dependent calcula- 
tion of the parameters fp, 5 2 and A are ambiguous in higher orders in a though 
the parametrization (6) itself is not. Another advantage of the proposed approach 
is a possibility to sum up many corrections with the help of a nonlocal operator 
expansion/ 1 zo 
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